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ABSTRACT 

We have carried out hydro- dynamical simulations to investigate the formation and 
evolution of protostar and circumstellar disks from the prestellar cloud. As the ini- 
tial state, we adopt the molecular cloud core with two non-dimensional parameters 
representing the thermal and rotational energies. With these parameters, we make 17 
models and calculate the cloud evolution ~ 10 4 years after the protostar formation. 
We find that early evolution of the star-disk system can be qualitatively classified into 
four modes: the massive disk, early fragmentation, late fragmentation, and protostar 
dominant modes. In the 'massive disk mode' to which the majority of models belong, 
the disk mass is greater than the protostellar mass for over 10 4 years and no frag- 
mentation occurs in the circumstellar disk. The collapsing cloud shows fragmentation 
before the protostar formation in the 'early fragmentation mode'. The circumstellar 
disk shows fragmentation after the protostar formation in the 'late fragmentation 
mode', in which the secondary star substantially gains its mass from the circumstellar 
disk after fragmentation and it has a mass comparable to that of the primary star. 
The protostellar mass rapidly increases and exceeds the circumstellar disk mass in the 
'protostar dominant mode'. This mode appears only when the initial molecular cloud 
core has a considerably small rotational energy. Comparison of our results with obser- 
vations indicates that a majority of protostar has a fairly massive disk during the main 
accretion phase: the circumstellar disk mass is comparable to or more massive than 
the protostar. It is expected that such a massive disk promotes gas-giant formation 
by gravitational instability in a subsequent evolution stage. 

Key words: star formation - circumstellar disk - methods: hydrodynamics - 
smoothed particle hydrodynamics - protoplanetary disk - planet formation 



1 INTRODUCTION 

There are two major mechanisms for gas-giant planet for- 
mation: one is the core accretion mechanism in which a mas- 
sive solid core forms first and the disk gas accretes onto the 
core (|Safronovlll969l : iGoldreich fc Wardlll973l : IPollack et al. 
1996), and the other is the gravitational instability (GI) 
mechanism in which the circumstell ar disk directly frag- 
ments into gas-giant planets via GI (|C amero ni l 19781 ). Re- 
cent discovery of extra-solar planets at a great distance from 
the central star such as HR 8799b, c, d and e (iMarois et al7\ 
120081 . l201Ch and GJ579b (|Thalmann et al. \\20od ) creates a 
new problem for the planet formation. It is difficult to form 
planets in the regions far from central stars according to the 
core accretion mechanism, because massive solid core for- 
mation before the dissipation of gaseo us disk seems to be 
difficult (jDodson- Robinson et al. \\200$ ). 



The gravitationally instability mechanism may be more 
plausible for the formation of these planets. Many stud- 
ies of disk fragmenta tion have be en done using either 
an analytic approach (Raflkov 2005) or numerical simula- 



tions (e.g., Stamatel los fc Whitworthll2007l ; ICai et al. 
"* "20061: 



Bolev et ' qZ.1 [2 006: Mei ia et oZ.1 120051 : [Pickett et al. 



2008; 



2003; 



Mem &; Batdl2010l V These efforts, however, seem to build 
a consensus that the planet formation by GI within ~ 50 
AU is highly difficult when the disk-to-stellar mass ratio is 
Mdisk/M s tar S, 0-1 whose ratio i s suggested by observations 
(see, e.g.. iKitamura et al. Il2002h . 



On the other hand, linutsuka et al. I (|201Qh showed that 
the circumstellar disk is comparable to or more massive 
than the protostar (Mdi s k/M st ar > 1) during the (early) 
main accretion phase (i.e., Class or Class I stages) and 
is highly gravitationally unstable. Recently, such massive 
disks were observed around very young protostars (e.g., 
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lEisner fe Carpenter! l200d : lEnoch et all [2009). Fragmenta- 
tion in such massive disks may potentially account for the 
formation of gas-giant planets with a wide separation such 
as HR 8799 or GJ579. Thus, we may have to focus on frag- 
mentation and gas-giant planet formation in a massive disk 
during the main accretion phase. There are several unsolved 
problems: how long such a massive disk exists, or whether 
the massive disk can fragment to form sub-stellar objects. 

So far, a few authors have investigated the early evolu- 
tion of the circumstellar disk and its fragmentation from 
prestellar core stage. Us ing the thin disk approximation, 
IVorobvov &; Bas"ul (|2010ah calculated the formation and evo- 
lution of the circumstellar disk in the molecular cloud core 
with a barotropic e quation of state (or incl uding cooling 
and heating effects, IVorobvov fc Basull2010bh . They found 
that a gravitationally unstable disk frequently appears in 
the star formation process and tends to show fragmentation 
during t he main accretion ph ase. With isothermal equation 
of state, iKratter et al. (2010) investigated the formation of 
the circumstellar disk in their three dimensional calculation. 
They also showed massive disk formation during the main 
accretion phase and concluded that the disk fragments if the 
disk-to- stellar mass ratio is greater than Mdisk/M sta r > 1. 
Using an adi a batic equation of state and radiative cooling, 
Walch et al. I (|2009h studied the circumstellar disk forma- 
tion with a somewhat coarse spatial resolution. They showed 
that a massive disk becomes hot during the early stage of 
the main accretion phase and fragmentation is suppressed. 

Machida et al. 1 (|2010h also investigat ed the di sk forma- 
tion with a finer spatial resolution than I Walch et a/.l ([2009) 
using a barotropic equation of state. In the collapsing cloud, 
the adiabatic heating dominates the radiative cooling at 
n rsj 10 10 cm -3 and the first (adiabatic) core with a size of 
^1AU forms before the pro t ostar formation (Larson 2001; 
iMasunaga &; InutsukafeoiOh . iMachida et al. I (|2010h showed 
that the first core formed before the protostar formation 
directly evolves into the circumstellar disk after the proto- 
star formation. Thus, to investigate the (early) formation of 
the disk, we need to resolve the first core that has a size 
of rsj 1 AU. In their calculation, however, they implicitly as- 
sumed the point-symmetric structure because they fixed the 
proto stellar location to the cent er of the computational do- 
main. iMichael &; Durisenl (|2010h pointed out that the stellar 
motion weakens GI activity in some degree, when the disk- 
to- stellar ratio is Mdi S k/M s tar ~ 0.1. Thus, a more adequate 
treatment may be necessary for the central protostar. 

All previous studies of circumstellar disk formation have 
shown that the circumstellar disk can become more mas- 
sive than a protostar during the mass accretion phase. How- 
ever, in such studies, the disk evolution was i nvestigated in a 
limit e d parameter range (iWalch et al. 2009; Machida et al. 



2010; IMachida fe M atsumoto 20lT|). Thus, we cannot con- 
clude whether such a massive disk that is a favorable site 
for gas-giant planet formation generally appears in the star 
formation process. To determine this, we need to investigate 
the disk formation resolving ^1AU structure in a wider pa- 
rameter range. 

The initial cloud for star formation is conventionally 
characterized by two parameters: the ratio of thermal a and 
rotational /3 energy toward the gravitational energy (for de- 
tailed description, see §2.2). With these parameters, many 
studies have calculated and classified the cloud evolution in 



the isother mal po ss 1993; Mivama 1984; Tsuribe Il2002h and 
adiab atic (|Cha &; Whitworthl 120031 : iMatsumoto &; Hanawal 
2003) gas contraction phases to investigate fragmentation 
before the protostar formation. However, the evolution and 
fragmentation of the circumstellar disk that is formed after 
the protostar formation has not been investigated with these 
parameters. 

In this study, with parameters of a and f3, we simu- 
late the formation and evolution of circumstellar disk from 
molecular cloud core until ~ 10 4 years after the protostar 
formation using a smoothed particle hydrodynamics code 
with sufficient spatial resolution. Parameters a and f3 signif- 
icantly influence the disk formation, because the accretion 
rate onto the circumstellar disk and the size and mass of the 
circumstellar disk are determined by these parameters (see, 
S}. We calculate the cloud evolution for 17 models using 
the barotropic equation of state. 

In 321 we describe the numerical method and initial con- 
ditions. The evolution of circumstellar disk with is presented 
in 33 We discuss our results and compare them with previ- 
ous works in 3H Finally, we summarize our results in §5\ 



2 NUMERICAL METHOD AND INITIAL 
CONDITION 

Our simulations are carried out using an SPH code newly 
developed for this work. The code includes an individual 
time-steps technique and uses Barnes-Hut tree algorithm to 
calculate self-gravity with opening angle — 0.5. We use 
adapt ive softening length according to IPrice &; Monaghanl 
(|2QQ7h . Arti ficial viscos i ty is i ncluded according to the pre- 
scriptio n bvlMonaghanl (|l997h with a v = 1 plus the Balsara 
switch (|Balsaralll995h . Our code is parallelized with MPI. 

To mimic the thermal evolution of molecular cloud, we 
use the barotropic equation of state as 



C 2 s,0P 



1+1 -e- 



2/5' 



(i) 



where c S) o = 190m s _1 and p c — 4xl0 _14 g cm -3 is adopted. 
As the initial state, we take a uniform density sphere with 
axisymmetric density perturbation. The density perturba- 
tion 5p is given as 5p = 0.01 x cos 20. The initial cloud is 
rigidly rotating. We parameterized the ratio of thermal and 
rotational energy to the gravitational energy of the initial 
cloud. They are characterized by two non-dimensional pa- 
rameters 



and 



^thermal _ 5Ro c t,0 
|£ ff ravl ~ 2GM ' 



3GM' 



(2) 



(3) 



where i?o, and M are the initial cloud radius, angu- 
lar velocity and mass of cloud core, respectively. We vary 
parameters a and /3 in the range of 0.4 < a < 0.8 and 
3 x 10~ 4 < P < 7 x 10 -2 . In this study, total cloud mass 
is fixed to 1M© in all models to compare the evolution of 
cloud with the same mass. In addition, we fixed the initial 
cloud temperature to 10 K. Thus, parameter a and initial 
cloud density are determined when the initial cloud radius 
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is given. Note that such treatment changes the gravitational 
energy of the cloud, After the cloud radius is fixed with arbi- 
trary a, parameter f3 is determined when the initial angular 
velocity Qo is given. The model name, parameters a and /3, 
initial cloud radius, initial angular velocity and initial cloud 
density are listed in Table 1. The initial cloud is modeled 
with about 500 000 SPH particles. 

To calculate the disk formation several 10 3 — 10 4 years 
after the protostar formation, we adopt a sink partic l e tech - 
nique according to the prescription by Bat e et all {l995). 
Starting from the prestellar core stage, the cloud evolution 
is calculated without sink particle. Then, we assume the pro- 
tostar formation and dynamically introduce a sink particle 
when the gas particle density exceeds the threshold den- 
sity, />sink = 4x 10 _9 g cm -3 . The threshold density roughly 
corresponds to th e densit y at w hich the second collapse be- 
gins. A s shown in [Larson] (|200lh and lMasunaga &; Inutsukal 
(2010|), the second collapse begins when the gas density 
reaches p ~ 10 _9 g cm -3 ~ /Wk at which the gas tem- 
perature exceeds T > 2 x 10 3 K and molecular hydrogen 
begins to dissociate. The protostar forms immediately after 
the second collapse. Thus, in this paper, we safely define the 
protostar formation epoch as that at which the gas density 
exceeds p > p s ink- 

To treat the gas accretion onto the sink particle after 
the creation of the sink particle, we set the accretion radius 
r aC c = 0.85AU. Then, within the accretion radius, we allow 
the gas accretion onto the sink particle when the following 
condition is fulfilled: (1) the gas particle density exceeds the 
accretion density p acc = 4 x 10 _11 g cm -3 , (2) it is grav- 
itationally bound to the sink particle and (3) its specific 
angular momentum is less than that required for it to form 
a circular orbit at r acc . We did not implement a boundary 
condition for sink particle. 



3 RESULTS 

As listed in Table 1, we calculated the cloud evolution for 
17 models in total, and classified them into four modes as 
follows: 

(i) Massive disk mode: the disk mass dominates the proto- 
stellar mass for over 10 4 years after the protostar formation. 

(ii) Early fragmentation mode: fragmentation occurs in 
the collapsing cloud before the protostar formation. 

(iii) Late fragmentation mode: fragmentation occurs in 
the circumstellar disk after the protostar formation. 

(iv) Protostar Dominant mode: the protostellar mass 
rapidly dominates the disk mass within 10 4 years after the 
protostar formation. 

The calculation results and its classification are summarized 
in Figure [TJ 

3.1 Massive Disk Mode 

Massive disk mode is indicated by the circles in Figure [TJ 
Figures [2] and [3] show the time evolution of the center of 
the cloud for model 8 (a = 0.6 and f3 = 7 x 10~ 3 ) that is 
a typical model for massive disk mode. The upper left and 
upper middle panels in these figures show snapshots before 
the protostar formation. A strong bar-like structure appears 



in the central high-density gas region in the top left panel in 
Figure [21 Then, the bar structure effectively transfers the an- 
gular momentum outward and the high-density gas region 
shrinks. As a result, a bimodal structure composed of the 
central high-density core and its surrounding disk appears 
as seen in the top middle panel. These figures clearly show 
the disk formation before the protostar formation. The top 
right panel shows the snapshot just when the second collapse 
occurs (i.e., the protostar forms). The time when the proto- 
star forms is 8.6 x 10 5 years from beginning. This roughly 
corresponds to the free-fall timescale of the initial cloud. By 
this epoch, the spiral arm has developed around the central 
object. The bottom left, middle, and right panels show the 
snapshots 1.4 x 10 2 , 1.0 x 10 3 and 1.0 x 10 4 years after the 
protostar formation. These figures show that the disk grad- 
ually increases its size with time retaining the global spiral 
structure. 

Figure U shows the mass evolution of the disk and pro- 
tostar. For this model, the disk mass is greater than the 
protostellar mass for more than 10 4 years after the proto- 
star formation. Such a disk is expected to be gravitationally 
unstable. To investigate the disk stability, the contours of 
Toomre's Q parameter at 1.4 x 10 2 and 1.0 x 10 3 years after 
the protostar formation (the same epoch of the bottom left 
and middle panels of Fig. [2} are plotted in Figure [5] The 
Toomre's Q parameter is described as 



where c s , k and E are the sound velocity and epicyclic fre- 
quency and surface density of the disk, respectively. For 
this model, the circumstellar disk did not show fragmen- 
tation even though the circumstellar disk has the region of 
Q < 1 at 1.4 x 10 2 years after the protostar formation. This 
is because GI cannot grow sufficiently fast in the circumstel- 
lar disk. The characteristic timescale of Toomre's analysis, 
r = 2c s /[GT,(l — Q 2 ) 1 ^ 2 ], is comparable to the orbital pe- 
riod of the disk (both are several thousand years). In this 
case, there are non-linear stabilize mechanisms against GI. 
As seen in the right panel of Figure [5] or bottom middle 
panel of Figure [2] the spiral arm globally redistribute the 
mass and angular momentum in a short duration, ~ 10 3 
years (< T orhi t). 

Furthermore, the disk- star configuration also dynami- 
cally changes in a short duration. Figure [6] shows the trajec- 
tory of the protostar for 8.0 x 10 3 years after the protostar 
formation. The asterisk indicates the position where the pro- 
tostar forms. The central protostar drifts toward the dense 
part of the spiral arm due to the gravitational interaction 
between the protostar and spiral arms. Since the strength of 
Keplerian shear is proportional to J^(rQ 2 ) oc r -3 , its radial 
dependency is very strong and the drift motion increases 
the shear stress of the dense region and suppresses frag- 
menting of the disk. As described in iGammid ([2001) , the 
thermal evolution of the disk also may play an important 
role to suppressing/promoting fragmentation. They showed 
that no fragmentation occurs when the cooling time is much 
longer than the orbital period. In our simulation, since we 
adopted the barotropic equation of state (see, §2.2), the disk 
evolves adiabatically (practically all the time, the disk mid- 
plane density is greater than p c — 4 x 10 _14 g cm -3 ). With 
these effects, the disk is stabilized as seen in the right panel 
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Table 1. Model parameters and Calculation Results. 
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Figure 1. The classification of simulation results on a- (3 plane. The circle, open triangle, filled triangle, cross indicate massive disk mode, 
early fragmentation mode, late fragmentation mode and protostar dominant mode, respectively. M and r ceri t are the mass accretion rate 
and centrifugal radius, respectively. 
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Figure 2. Time sequence of the logarithm of the face-on surface density before and after the protostar formation for model 8 (a = 0.6 and 
P = 7 x 10 -3 ). Top left and top middle panels show the snapshot about 1.1 x 10 3 and 6.9 x 10 2 years before the protostar formation. Top 
right panel shows the snapshot just when the protostar forms. The bottom left, middle and right show the snapshots 1.4 x 10 2 , 1.0 x 10 3 
and 1.0 x 10 4 years after the protostar formation. The elapsed time from beginning is described in each panel. 
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Figure 4. Mass evolution of the protostar (dashed line) and the 
disk (solid line) for model 8. 



of figure [5] For disk fragmentation, the mass in- fall must be 
sufficiently fast to overcome these non-linear back reactions. 



3.2 Early Fragmentation Mode 

The open triangles in Figure [1] indicate " Early fragmenta- 
tion mode," in which the collapsing cloud undergoes frag- 
mentation before the protostar formation. Figure [8] shows 
the time evolution of the center of the cloud for model 11 
(a = 0.4 and f3 = 7 x 10 -2 ). The figure shows that the ring 
pattern arises (top left panel) and it fragments to form four 
clumps (or protostars). These protostars have roughly the 
same mass when they form. After fragmentation, they inter- 
act with each other. Finally, a quadruple stellar system com- 
posed of two close binary appears as seen in the lower right 
panel of Figure [8] The conditio n of the ring like structur e 
was substantially investigated in ICha &; Whitworthl (|2QQ3h . 
Since they fixed parameter a (but changed the rotational 
low), we cannot quantitatively compare their results with 
ours. However, our result is not contradict with them. 

In models showing fragmentation before the protostar 
formation (i.e., early fragmentation mode), fragmentation 
pattern qualitatively changes as the rotation energy of the 
initial cloud decreases. Figure [9] shows the time evolution 
of the center of the cloud for the model 12 (a = 0.4 and 
P = 1.0 X 10 -2 ). A clear bar- like structure arises before the 
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Figure 5. The contours of Toomre's Q parameter at 1.4 x 10 2 years (left) and 1.0 x 10 3 years (right) after the protostar formation. 




Figure 6. The protostellar trajectory for 8.0 x 10 3 years after the 
protostar formation. The asterisk indicates the position where the 
protostar forms. 



protostar formation for model 12. The bar fragments into 
two clumps. Then, a binary system appears around the cen- 
ter of the collapsing cloud (upper middle panel of Fig. [9j. 
Furthermore, the remnant of the bar fragments and form 
two extra protostars about 10 3 years after the first fragmen- 
tation, as seen in the upper right panel of Figure [9] Then, 
one protostar is ejected from the central region by gravi- 
tational interaction among protostars, and a triple stellar 
system remains around the center of the cloud at the end of 
the calculation. 

Fragmentation and subsequent evolution of fragments 
are very complicated in models belonging to the early 
fragmentation mode. However, since fragmentation before 
the protostar fo rmation was well investigated in p revi- 
ous studies (see, Tsuribe 2002, iGoodwin et al I l2007l and 




Bodenhei mer et al. l bOOOl for review), we do not comment it 



any more in this paper. 



1000 2000 3000 4000 5000 6000 7000 8000 9000 
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Figure 7. Mass evolution of the primary star (solid line) and the 
secondary star (dashed line). 



3.3 Late Fragmentation Mode 

The filled triangles in Figure [1] indicates the "Late fragmen- 
tation mode," in which fragmentation occurs in the circum- 
stellar disk after the protostar formation. Models 13 and 
14 belong to this mode. Figure [10] shows the snapshots for 
model 14 (a = 0.4 and /3 = 7 x 10" 3 ). As seen in the mas- 
sive disk mode, after a strong bar develops (top left panel), 
bimodal structure composed of the central high-density gas 
region and disk-like structure appears even for model 14 
(top middle panel). However, fragmentation occurs in the 
disk before the spiral arm sufficiently develops (top right 
panel) in this model, while no fragmentation occurs in the 
models belonging to the massive disk mode (see, Figs. 2 
and 3). After fragmentation, the secondary object appears 
in the region ~ 30 AU far from the primary star (bottom 
left panel) ~ 6 x 10 2 years after the primary star formation. 
Thus, Figure [10] indicates that fragmentation can occur in 
the circumstellar disk during the main accretion phase if the 
mass in-fall is sufficiently fast (with smaller a, see §5$. 

The mass evolution of two protostars is plotted in Fig- 
ure [7] in which the secondary star increases its mass in a 
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way similar to that of the primary star. At the end of the 
calculation, the mass of secondary star is comparable to that 
of the primary star. Thus, this system evolves into a binary 
stellar system, not a star-planet system. 



3.4 Protostar Dominant Mode 

The crosses in Figure [1] indicate "protostar dominant mode", 
in which the protostar rapidly increases its mass and ex- 
ceeds the circumstellar disk mass within 10 4 years. For this 
mode, the protostellar mass exceeds the disk mass within 
10 4 years after the protostar formation. Figures [11] and IT2l 
show the time evolution of the center of the cloud for model 
16 (a = 0.4 and (3 = 10" 3 ) that is a typical model for the 
protostar dominant mode. The top left and middle panels 
show the snapshots before the protostar formation for model 
16. Unlike the massive disk mode, the first core keeps an 
almost axisymmetric structure without developing the non- 
axisymmetric perturbation. The second collapse occurs at 
4.4 x 10 4 years from the beginning. Even after the protostar 
formation, the axisymmetric structure is maintained until 
the disk radius sufficiently grows (r > 10 AU). 

To investigate the disk stability for model 11, the 
Toomre's Q parameter at 1.1 x 10 3 years after the proto- 
star formation (the same epoch as the bottom middle of 
Fig. fTT]) is plotted in Figure [T3l At this epoch, the disk mass 
is still greater than the protostellar mass (see, Fig. ll4p . Nev- 
ertheless, the Q parameter is greater than unity in the whole 
region of the disk due to compactness of the disk. 

The mass of the protostar and circumstellar disk is plot- 
ted against time after the protostar formation in Figure IT4l 
The figure indicates that the accretion rate onto the pro- 
tostar is higher than that of the circumstellar disk, and 
the protostellar mass exceeds the circumstellar disk mass 
at rsj 2.0 x 10 3 years after the protostar formation. 



circumstellar disk is typically about 10 4 years. Thus, the ra- 
diative cooling may play an important role in investigating 
the thermal evolution of the disk and its fragmentation over 
10 4 years. In this study, however, we showed that the evolu- 
tion of the circumstellar disk can be qualitatively classified 
into four modes during very early stages of the star forma- 
tion (t < 10 4 years). Thus, we expect that our classification 
does not change qualitatively even when the radiative effects 
are included. However, the radiative cooling of disk is im- 
portant to investigate the disk evolution and fragmentation 
for t > 10 4 years. 

4.2 Treatment of Protostar and Sink 

In this study, no model shows the formation of a star-planet 
system during the main accretion phase because the sec- 
ondary object continues to increase its mass and finally ex- 
ceeds the h ydrogen- burning limit (M > 0. 08 Mp). On the 
other hand, IVorobvov &; Basul (|2010al ) and iMachida et al. I 
(2010) showed the formation of a star-planet system during 
the main accretion phase. The difference is thought to be 

caused by treatment of the p rotosta r and sink. 

In IVorobvov &; Basul (|2010ah and IMachida et al. I 

(|2010l ) , the protostar (or sink cell) is fixed at the center of the 
computational domain. It is expected that such treatment 
promotes fragmentation. In reality, the density fluctuation 
arising around the protostar can cancel out by movement of 
the protostar. Thus, fragmentation tends to occur when the 
protostar is fixed. 

In addit i on, IVorobvov fe Basul (|2010aT > and 
Machida et al. (2010) did not impose the sink on fragments 
formed in the circumstellar disk. Instead, they suppressed 
further collapse of fragments with adiabatic equation 
of state. Such treatment decreases the mass accretion 
onto fragments in some degree. On the other hand, our 
sink treatment may overestimate the mass accretion onto 
fragments or protostar. 



4 DISCUSSION 



4.3 Comparison with Previous Works 



4.1 Effect of Radiative Cooling 

In this study, we used the barotropic equation of state 
(eq. [1]), in which the gas in the disk behaves adiabati- 
cally when the disk mid-plane density is greater than p c = 
4x 10 -14 cm -3 . However , in reality, the disk radiatively cools 
with time. According to iRafikovl (|2005r h we can roughly es- 
timate the cooling time of the disk t coo i as 

„ /to 



1 2crT 4 



1.2 x 10 4 



200g cm- 



(—)' 

V iook J 



2.3 



lOcnr 



where E, a, T, r and fl are the surface density, Stephan- 
Boltzmann constant, mid-plane temperature, optical depth 
of the disk and mean molecular weight, respectively. In equa- 
tion ([5]), we estimate the vertical optical depth f(r) of the 
disk as f(r) — (1 + t 2 )/t. In addition, we assume the optical 
depth as r~ «E/2, where k is the opacity. We assume that 
the disk is optically thick, and derive RHS of eq J5]) with 
fir) ~ r. 

Equation J5]) indicates that the cooling timescale of the 



So far, few authors have investigated the evolution of cir- 
cumstellar disk from m olecular cloud c o re. W ith an isother- 
mal equation of state, iKratter et al. showed that 
fragmentation occurs in the circumstellar disk with a wide 
parameter space during the early stage in the main accre- 
tion phase and claimed that fragmentation frequently occurs 
when the disk-to-stellar mass ratio is greater than unity. 
However, their assumption of isothermality seems not to 
be valid for the disk evolution of the early main accretion 
phase, because the gas becomes opaque and behaves adia- 
batically when the gas density exce eds the critical density 
-e^oj (yei((S^ 13 - 10 _14 g cm -3 (e.g., iMasunaga Inutsukal 
§Oi0). In addition, the radiative cooling can affect the disk 
evolution ~ 10 4 years after the protostar formation as de- 
scribed in £14.11 Since the adiabatic equation of state stabi- 
lizes the circumstellar disk, fragmentation barely occurs in 
our calculation eve n when the dis k -to-st ellar mass ratio ex- 
ceeds unity. Thus, Kra tter et al. may overestimate 
the fragmentation condition, while our calculation may un- 
dere stimate it because o f lack of radiative cooling. 

Walch et al. (2009) also studied the circumstellar disk 
formation with the adiabatic equation of state and radiative 
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Figure 8. Time sequence of the logarithm of the face-on surface density for model 11 (a 
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Figure 9. Same as Fig. El but for model 12 (a = 0.4 and (3 = 1.0 X 10~ 2 ). 



cooling. They showed no fragmentation in the circumstellar 
disk, because the disk becomes very hot during the early 
stage of main accretion phase. Their spatial resolution is, 
however, somewhat coarse; the minimum smoothing length 
of their study is /i m i n = 2 AU. On the other hand, /i m i n ~ 
0.3 AU in our simulations. In addition, they restricted their 
initial conditions to rapidly rotating cases (£ > 10~ 2 ) to 
investigate disk evolution before the central density becomes 
high. The observation suggested that molecular cloud cores 
have the rotational energy of 10~ 4 < 6 < 0.07 with a typical 
value of p ~ 0.02 (|Caselli et al Il2002h . Thus, they studied 
the cloud evolution in the limited parameter range. 



5 SUMMARY 

We have carried out hydro- dynamical simulation to investi- 
gate the evolution of the circumstellar disk with two non- 
dimensional parameters representing the thermal and rota- 
tional energy of the initial cloud. The thermal energy a is 
related to the mass accretion r ate onto the circumst ellar disk 
as Mdisk = a" 3/2 cfG _1 fsee. iMachida et al. Il201lh . Thus, 
smaller a provides a high accretion rate onto the circum- 
stellar disk and vice versa. On the other hand, the initial 
rotational energy that is represented by parameter f3 is re- 
lated to the disk radius. The centrifugal radius of the initial 
cloud is related to /3 as r ce nt = 3RoP, where Ro is the initial 
cloud radius. Thus, with larger f3, the cloud forms a larger 
disk in the main accretion phase. In other words, with larger 
/3, a large fraction of the in-falling matter accretes onto the 
disk, rather than directly onto the primary protostar. As a 
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Figure 11. Time sequence of the logarithm of the face-on surface density for model 16 (a ■ 
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Figure 12. Same as figure [TT] but edge-on snapshots are shown. 
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Figure 14. Mass evolution of the protostar (dashed line) and the 
disk (solid line) for model 16. 



result, smaller a and larger f3 increase disk surface density 
and makes a gravitationally unstable disk. 

On the other hand, the non-axisymmetric structure 
arose in such unstable disk can stabilize the disk, because 
it can redistribute the angular momentum and promote 
mass accretion onto the protostar. Thus, no fragmenta- 
tion occurs when a strong non-axisymmetricity grows and 
transfers sufficient angular momentum outward. By con- 
trast, the disk becomes highly gravitationally unstable and 
shows fragmentation when non-axisymmetric structure does 
not grow sufficiently or when the growth timescale of the 
non-axisymmetricity is much longer than the disk growth 
timescale. Thus, fragmentation condition depends also on 
the growth of the non-axisymmetriciy, which is closely re- 
lated to parameters a and f3 because they determine the 
evolution of the mass and angular momentum of the disk. 

With parameters a and j3, we found that the disk evo- 
lution is qualitatively classified into four modes: protostar 
dominant, massive disk, early fragmentation and late frag- 



mentation modes. The schematic classification of the cir- 
cumstellar disk is shown in Figure [1] that covers a wide pa- 
rameter range of the rotati onal energy suppo rted by obser- 
vations (1(T 4 < p < 0.07; ICaselli et aflbOO^ . We describe 
each mode in the following. 



• Protostar Dominant Mode: For protostar dominant 
mode, the protostellar mass exceeds the circumstellar disk 
mass within 10 4 years after protostar formation. This mode 
appears in the range of ft < (1 — 3) x 10 -3 when the protostar 
forms, the circumstellar disk is more massive than the pro- 
tostar. However, for this mode, since the mass increase rate 
of the protostar is larger than that of the circumstellar disk, 
the protostellar mass exceeds the circumstellar disk mass in 
~ 10 4 years after the protostar formation, as described in 

go 

• Massive Disk Mode: For the massive disk mode, the 
circumstellar disk is more massive than the protostar over 
10 4 years after protostar formation. As seen in Figure [T] 
this mode appears when the initial cloud has larger rota- 
tional and thermal energies. In our calculation, a majority 
of models (8 out of 17 models) belongs to this mode. This in- 
dicates that the circumstellar disk that is more massive than 
the protostar frequently appears in the star formation pro- 
cess. The circumstellar disk is self-regulated for this mode. 
When the circumstellar disk becomes massive, it decreases 
its mass because the non-axisymmetric structure develops 
and promotes the mass accretion onto the protostar. On the 
other hand, when the disk is relatively less massive, the disk 
mass increases because a relatively stable disk acquires its 
mass from the in-falling envelope without effective angular 
momentum transfer. 

• Early Fragmentation Mode: For early fragmenta- 
tion mode, fragmentation occurs before the protostar for- 
mation. As seen in Figure [T] this mode appears when the 
initial cloud has large rotational but small thermal energies. 
If the initial rotational energy is large enough, the gas can- 
not condense to the center, and a ring like structure appears 
as seen in figure [8] This structure is unstable against per- 
turbation and it fragments into protostars. As the rotational 
energy decreases, the gas can condense to the center and bar 
like structure develops. However, unlike in the massive disk 
mode, since the bar cannot transfer angular momentum out- 
ward fast enough, fragmentation occurs and binary system 
appears as seen in Figure [9] 

• Late Fragmentation Mode: For late fragmentation 
mode, fragmentation occurs after the protostar formation. 
The rotational energy for this mode is smaller than that 
for early fragmentation mode. Thus, without fragmentation 
before the protostar formation, the central bar like structure 
shrinks transferring angular momentum outward and form 
a single protostar and a disk. 



In this study, since we did not include proper radiation 
treatment, we cannot discuss whether a gas-giant planet can 
form in such a massive disk. However, such a massive disk is 
a plausible site for the gas-giant formation by GI. To deter- 
mine whether a gas-giant planet forms in the circumstellar 
disk by GI, we need to perform more long-term calculation 
including radiative effects in the future. 
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